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ABSTRACT 

We have performed detailed numerical calculations of the non-conservative evolution of close X-ray binary 
systems with intermediate-mass (2.0 - 6.OM0) donor stars and a \3Mq accreting neutron star We calculated 
the thermal response of the donor star to mass loss, in order to determine its stability and follow the evolution 
of the mass transfer. Under the assumption of the "isotropic re-emission model" we demonstrate that in many 
cases it is possible for the binary to prevent a spiral-in and survive a highly super-Eddington mass-transfer phase 
(1 ^ M /Mem < 10'') on a sub-thermal timescale, if the convective envelope of the donor star is not too deep. 
These systems thus provide a new formation channel for binary millisecond pulsars with heavy CO white dwarfs 
and relatively short orbital periods (3-50 days). However, we conclude that to produce a binary pulsar with a 
O-Ne-Mg white dwarf or PqA ^ 1 day (e.g. PSR B0655H-64) the above scenario does not work, and a spiral-in 
phase is still considered the most plausible scenario for the formation of such a system. 

Subject headings: stars: evolution, mass-loss, neutron — white dwarfs — binaries: close 
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1. INTRODUCTION 

Recently a large number of binary millisecond pulsars (BM- 
SPs) with relatively heavy white dwarf (WD) companions have 
been reported. These pulsars form a distinct class of BM- 
SPs (cf. Table 1) which are characterized by relatively slow 
spin periods (fspin — 10-200 ms) and high period derivatives: 
10"^° < Pspin < 10"'^. It has been suggested that such BMSPs 
evolved through a common envelope and spiral-in phase (e.g. 
van den Heuvel 1994). This gives a natural explanation for 
their close orbits and the presence of a relatively heavy CO/O- 
Ne-Mg WD, if the mass transfer was initiated while the donor 
star (the progenitor of the WD) ascended the AGB. However, 
it has been argued that a neutron star engulfed in a common 
envelope might experience hypercritical accretion and thereby 
collapse into a black hole (e.g. Chevalier 1993; Brown 1995). 
If this picture is correct then these BMSPs can not have formed 
in a common envelope and spiral-in phase. 

Here we investigate an alternative scenario for producing the 
mildly recycled BMSPs with He or CO WD in close orbits, in 
which a 2-6 Mq donor star, with a non- (or partly) convective 
envelope, transfers mass on a sub-thermal timescale and yet in 
a dynamically stable mode. 

2. STABILITY CRITERIA AND MODE OF MASS TRANSFER 

The stability and nature of the mass transfer is very impor- 
tant in binary stellar evolution. It depends on the response of the 
mass-losing donor star and of the Roche-lobe (e.g. Paczynski 
1976; Soberman, Phinney & van den Heuvel 1997). The mass 
transfer is stable as long as the donor star's Roche-lobe contin- 
ues to enclose the star. Otherwise it is unstable and proceeds on 
the shortest unstable timescale. 

As long as the mass of the donor, M2 is less than I.8M0 
the mass transfer will be dynamically stable for all initial or- 
bital periods (e.g. Tauris & Savonije 1999). These LMXBs 
are the progenitors of the BMSPs with a helium WD compan- 
ion. The observational absence of X-ray binaries with Roche- 
lobe filling companions more massive than ^ 2Mq has been at- 

' Center for High-Energy Astrophysics, and Astronomical Institute "Anton 
Netherlands; tauris@astro.uva.nl, edvdh@astro.uva.nl, gertjan@astro.uva.nl 



tributed to their inability to transfer mass in such a stable mode 
that the system becomes a persistent long-lived X-ray source 
(van den Heuvel 1975; Kalogera & Webbink 1996). Below we 
investigate for these systems how the stability of the RLO de- 
pends on the evolutionary status of the donor (and hence the 
orbital period) at the onset of mass transfer 

2.1. Numerical computations 

We have calculated the evolution of a large number of X-ray 
binaries with a donor star of mass 2 < M2 /Mq < 6 and a 1 .3Mq 
accreting neutron star Both the radius of the donor star as well 
as its Roche-lobe are functions of time and mass (as a conse- 
quence of nuclear burning, magnetic braking and other tidal 
spin-orbit couplings). We used an updated version of Eggle- 
ton's numerical computer code (Pols et al. 1998) to keep track 
of the stellar evolution and included a number of binary inter- 
actions to carefully follow the details of the mass-transfer pro- 
cess. For all donor stars considered here we assumed a chem- 
ical composition of X=0.70 and Z=0.02, and a mixing-length 
parameter of a = I /Hp = 2.0. We refer to Tauris & Savonije 
(1999) for a detailed description of our computer code. 

2.2. Highly super-Eddington mass transfer 

The maximum accretion rate onto a neutron star is given ap- 
proximately by the Eddington limit for spherical accretion of a 
hydrogen gas, MEdd = 1.5 x 10~^ Mq yr"'. If the mass-transfer 
rate from the donor star, M2 is larger than this limit, radiation 
pressure from the accreted material will cause the infalling mat- 
ter to be ejected from the system at a rate: 

|M I = IM2 1 - MEdd ^ \M2 1 if M2 > Mfidd 
In systems with very large mass-transfer rates, matter piles 
up around the neutron star and presumably forms a growing, 
bloated cloud engulfing a large fraction of the accretion disk. A 
system will only avoid a spiral-in if it manages to evaporate 
the bulk of the transferred matter via the liberated accretion 
energy. This would require the radius of the accretion cloud, 
Tci > ^Ns|M2|/MEdd in order for the liberated accretion en- 
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ergy to eject the transferee! material (~ O.lMEddC^ ^ 5 A/2 Vg^g, 
where v^^^ = 2GMns/''ci; ^ns is the radius of the neutron star). 
If the material which is to be ejected comes closer to the neu- 
tron star it will have too much negative binding energy in order 
for the liberated accretion energy to expel it. At the same time 
fci must be smaller than the Roche-lobe radius of the neutron 
star during the entire evolution, if formation of a common en- 
velope (CE) is to be avoided^. 

A simple isotropic re-emission model will approximately re- 
main valid for our scenario. In this model it is assumed that 
matter flows over conservatively from the donor star to the 
vicinity of the neutron star before it is ejected with the spe- 
cific orbital angular momentum of the neutron star. Assuming 
this to be the case we find that, even for extremely high mass- 
transfer rates (IM2I > lO^MEdd). the system can avoid a CE and 
spiral-in evolution. 

3. RESULTS 

3.1. A case study: M2 = 4.0Mq and Pa,-b = 4.0 days 

In Fig. 1 we show the evolution of a binary initially consist- 
ing of a neutron star and a zero age main-sequence companion 
star with masses Mns = 1.3Mq and M2 = 4.0Mq, respectively, 
and initial orbital period PqA = 4.0 days. 

At the age of f = 176.6 Myr the companion has evolved to 
fill its Roche-lobe (Ri = 8.95/?©; Jeff = 9550 K) and rapid mass 
transfer is initiated (A). The donor star has just evolved past the 
MS hook in the HR-diagram and is burning hydrogen in a shell 
around a O.56M0 hehum core. Prior to the mass-transfer phase, 
a radiation-driven wind (IM2I ~ 4 x 10"'" Mq yr"') has caused 
the donor to decrease its mass slightly (M2 = 3.99Mq) and con- 
sequently resulted in a slight widening of the orbit (PqA = 4.02 
days). Once the donor fills its Roche-lobe it is seen to lose mass 
at a very high rate of |M2 1 ^ 4 x 10"^ Mq yr"' =2.7x10^ Mudd- 
At this stage the donor has only developed a very thin convec- 
tive envelope of size Zcnv = O.OISRq, so its envelope is still 
radiative and will therefore shrink in response to mass loss. At 
f = 176.7 Myr its radius has decreased to a minimum value of 
3.38/?0, but now Z^om = 0.97 /J©. At this point (s) the donor 
has a mass of 1 .76M0, TeB = 5640 K, and Porb = 1 -59 days. The 
donor expands again, but shortly thereafter, its rate of expansion 
slows down causing IM2I to decrease to ^ lOMndd- 

The mass transfer ceases (B) when the donor has an age of 

178.3 Myr. At fliis stage Porb = 8.11 days, R = 6.68/?©, Zcom = 
0.07 /?0 and Teff = 12700 K. The mass of the donor is O.6I8M0. 
It still has a O.56M0 helium core, but now only a O.O6M0 en- 
velope consisting of 16% H and 82% He. The mass-transfer 
phase (A-B) lasts relatively short: ?x = 1 -7 Myr, and hence the 
neutron star will only accrete: AMns = fx Afedd = 0.03 M©. This 
leads to relative large values of Pspin and Pspin for the (mildly) re- 
cycled pulsar. As the mass-transfer rate is always highly super- 
Eddington during the RLO, and the accreting neutron star wiU 
be enshrouded by a thick (bloated) disk, it is doubtful whether it 
will be observable as an X-ray binary during this phase - except 
very briefly just at the onset and near the end of the RLO. 

We followed the evolution of the donor star further on. The 
donor continues to burn hydrogen in its light envelope. At / = 

186.4 Myr (f) the helium burning is finally ignited (Lue/Ln > 
10) in the core which now has a mass of O.596M0. After 70 
Myr (f = 253.8 Myr since the ZAMS) the core-hehum burning 



is exhausted (g). The O.6O2M0 core then has a chemical com- 
position of 19% C, 79% O and 2% Ne. It is surrounded by a 
O.OI6M0 envelope (16% H, 82% He, 1% N14). The central 
density is pc = 4.13 x 10^ g cm"^ and R2 = O.IIRq. From here 
on the star contracts and settles as a hot CO white dwarf. 
We have now demonstrated a scenario for producing a BMSP 
with the same characteristics as those listed in Table 1. 

3.2. The Porb-MwD diagram 

In Fig. 2 we have plotted the calculated final orbital periods 
as a function of the mass of the white dwarf companion (the 
remnant of the donor) for a given initial mass of the donor star. 
The values for the BMSPs which originated from a binary with 
a low-mass companion (the former LMXBs with M2 < 1 .SM© 
located on the upper branch) are taken from Tauris & Savonije 
(1999). These white dwarfs are expected to be helium WD - 
unless the initial orbital period was very large (PqA > 150 days) 
so a relatively heavy helium core developed prior to the RLO, 
in which case the helium core later ignited forming a CO WD. 
The final product of X-ray binaries with M2 > 2Mq are seen to 
deviate significantly from the low-mass branch. The reason is 
the former systems had a large mass ratio, q = M2/MNS which 
caused the binary separation to shrink initially upon mass trans- 
fer- cf. Fig. 1. Such systems only "survive" the mass-transfer 
phase if the envelope of the donor is radiative or slightly con- 
vective. This sets an upper limit on the initial orbital period for 
a given system. If the donor is in a wide binary it develops a 
deep convective envelope prior to filhng its Roche-lobe and it 
will therefore expand rapidly in response to mass loss which, in 
combination with the orbital shrinking, will result in the forma- 
tion of a CE and a (tidally unstable) spiral-in evolution. 
In Fig. 3 we show how the final orbital period and the mass of 
the WD depends on the initial orbital period for a binary with 
M2 = 4.OM0. We notice that the question of initiating RLO be- 
fore or after the termination of hydrogen core burning (case A 
or early case B, respectively) is important for these relations. 
For initial forb < 2.4 days (case A RLO), P^b decreases with 
increasing Porb- The reason is simply that in these systems the 
donor star is still on the the main-sequence and the mass of 
its helium core, at the onset of RLO, increases strongly with 
forb and therefore the amount of material to be transferred (the 
donor's envelope) decreases with Porb- Since the orbit widens 
efficiently near the end of the mass transfer, when the mass ra- 
tio between donor and accretor has been inverted (cf. Fig. 1), 
PqjI, will also decrease as a function of initial PqA- However, for 
Porb > 2.4 days (early case B RLO) the final orbital period in- 
creases with initial orbital period as expected - the core mass of 
the donor only increases slightly (due to hydrogen shell burn- 
ing) as a function of initial Porb- 

3.3. The initial (M2,Poib) parameter space 

In Fig. 4 we outline the results of our work in a diagram 
showing the fate of a binary as a function of its initial Porb 
and the value of M2. We conclude that X-ray binaries with 
2 ^ M2/M0 < 6 can avoid a spiral-in and CE evolution if Porb 
is between 1-20 days, depending on M2. If the initial Pqa is too 
short, the systems will obviously enter a CE phase, since these 
system always decrease their orbital separation when the mass 

^Note, that if no efficient cooling processes are present in the accretion disk then the incoming matter retains its net (positive) energy and is easily ejected in the 
form of a wind from the disk (Narayan & Yi 1995; Blandford & Begelman 1999). Even if the arriving gas is able to cool, interactions between the released radiation 
from this process and the infalling gas may also help to eject the matter. In both cases can be smaller than estimated above. 
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transfer is initiated^. On the other hand, if the initial Porh is too 
large the donor develops a deep convective envelope prior to 
RLO and a runaway mass-transfer event is unavoidable, also 
leading to a CE formation. For systems with M2 <, 1.8Mq 
and initial foit < 1 day the outcome is a BMSP with an ultra 
low-mass degenerate hydrogen star (e.g. PSR J205 1-0827, cf. 
Ergma, Sarna & Antipova 1998). 

4. DISCUSSION 

We have now demonstrated how to form a BMSP with a 
relatively heavy (He or CO) WD companion without evolving 
through a CE phase. If a substantial fraction of BMSPs have 
evolved through a phase with super-Eddington mass transfer 
on a sub-thermal timescale (a few Myr), this will eliminate the 
need for a long X-ray phase. This would therefore help solving 
the birthrate problem between BMSPs and LMXBs (Kulkarni 
& Narayan 1988) for systems with P^^.^^ < 50 days. 
It has recently been suggested (Podsiadlowski & Rappaport 
1999; King & Ritter 1999) that Cygnus X-2 descended from 
an intermediate-mass X-ray binary via a scenario which resem- 
bles the one described here. We confirm that Cygnus X-2 is a 



progenitor candidate for a BMSP with a heavy WD. 

It is seen from Fig. 2 that we can not reproduce the systems 
with very massive O-Ne-Mg WD or the short orbital periods 
(<^ 3 days) observed in some systems with a CO WD. We there- 
fore conclude that these binaries most likely evolved through a 
CE phase where frictional torques were responsible for their 
present short P^^^, (cf. thin lines in Fig. 2 and gray area in 
Fig. 4). These systems therefore seem to originate from bina- 
ries which initially had a relatively large Porb and case C RLO - 
otherwise if Po^b was small the stellar components would have 
coalesced either in the spiral-in process, or as a result of gravita- 
tional wave radiation shortly thereafter (typically within 1 Gyr 
for systems surviving case B RLO and spiral-in). 

We thank the Parkes Multibeam Survey team and the Swin- 
burne Pulsar Group for releasing binary parameters prior to 
publication. We appreciate comments from Ene Ergma on the 
issue of forming BHWD systems. T.M.T. would like to thank 
Gerry Brown for discussions and hospitality at Stony Brook. 
T.M.T. acknowledges the receipt of a Marie Curie Research 
Grant from the European Commission. 
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'in these narrow binaries the amount of available orbital energy (a possible energy source for providing the outward ejection of the envelope) is small and hence 
the neutron star is most likely to spiral in toward the (unevolved) core of the donor, forming a Thome-Zykow-like object. In that case the neutron star will probably 
undergo hypercritical accretion and collapse into a black hole. 
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Table 1 

Observed Pulsars with a "heavy" WD Companion 



PSR 


Porb 

(days) 


f 

(Mo) 


(Mq) 


^spin 

(ms) 


^spin 




J1904+04* 


15.75 


0.0046 


0.27 


71.1 






J1810-2005* 


15.01 


0.0085 


0.34 


32.8 


1.3 X 10" 


-19 


J1453-58* 


12.42 


0.13 


1.07 


45.3 






J0621+1002 


8.319 


0.0271 


0.540 


28.9 


< 8 X 10- 


-20 


J1022+1001 


7.805 


0.0833 


0.872 


16.5 


4.2 X 10" 


-20 


J2145-0750 


6.839 


0.0242 


0.515 


16.1 


2.9 X 10" 


-20 


J1603-7202 


6.309 


0.00881 


0.346 


14.8 


1.4 X 10" 


-20 


J1157-5112** 


3.507 


0.2546 


>1.20 


43.6 


< 9 X 10' 


-19 


J1232-6501* 
J 1435-60* 
B0655+64 
J1756-5322** 


1.863 
1.355 
1.029 
0.453 


0.0014 
0.14 
0.0714 
0.0475 


0.175 
1.10 
0.814 
0.683 


88.3 
9.35 
196 
8.87 


1.0 X 10" 
6.9 X 10" 


-18 
-19 



*New pulsar, Parkes Multibeam Survey (Manchester et al. 1999). 
**New pulsar, Edwards et al. (1999). 

Note. — is calculated assuming Mns = 1.4M0 and i = 60° . 
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Fig. 1 . — The evolution of an X-ray binary with M2 = 4.0 Mq and /"orb = 4.0 days. The left panel shows the evolution of Porb as a function of M2 (time is increasing 
to the right). The central panel gives the mass-loss rate of the donor as a function of its age since the ZAMS. The right panel shows the evolution of the mass-losing 
donor (solid line) in an HR-diagram. The dotted line represents the evolutionary track of a single 4.QMq star. The letters in the different panels correspond to one 
another at a given evolutionary epoch - see text for further explanation. 

Fig. 2. — The final P(„b as a function of WD mass for different BMSPs. Next to each curve is given the initial mass of the donor star (the progenitor of the WD) 
used in our evolutionary calculations. The free parameter in each curve is the initial P(„b (at the onset of the RLO). The curves in gray color represent the formaion 
of BMSPs with helium WD, while the black curves are BMSPs with CO WD. The open circles on some of the curves indicate the transition from case A to early 
case B RLO mass transfer (i.e. whether or not the donor burned hydrogen in the core at the onset of the RLO, Kippenhahn & Weigert 1990). The thin lines show the 
calculated parameters for systems which evolved through a CE and spiral-in phase scenario assuming an efficiency parameter of rjc^X =1.0 (e.g. van den Heuvel 
1994). The 12 observed BMSPs with a "heavy" WD companion are marked with a star, see Table L 

Fig. 3. — The dependence of final orbital period (top) and mass of the WD (bottom), on the initial orbital period, Pg^t,. 

Fig. 4. — This plot illustrates the allowed parameter space (white area) for producing BMSPs without evolving through a CE phase. If M2 > I.SMq and the 
donor has a deep convective envelope at the onset of mass transfer (i.e. P„ib is large) the system will evolve into a CE and spiral-in phase. This is also the case if the 
initial period is very short and M2 > I.8M0. In the latter case the neutron star may collapse into a black hole. 
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